Radiative decays of pseudoscalar (P) and vector (V) mesons and the process e + e — > rfp 
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Radiative decays of pseudoscalar and vector mesons are calculated in the framework of the chiral 
Nambu-Jona-Lasinio (NJL) model. We use the amplitude for triangle quark loops of anomalous 
type. In evaluating these loop integrals we use two methods. In the first one, we neglect the 
dependence of external momenta by reproducing the Wess-Zumino-Witten terms of effective chiral 
meson Lagrangian. In the second method, we take into account the momentum dependence of loop 
integrals omitting their imaginary part. This makes it possible to allow for quark confinement. As 
applied both the methods is in qualitative agreement with each other and with experimental data. 
The second method allows us to describe the electron-positron annihilation with production of rj 
and p mesons in the center of mass energy range from 1.6 to 3.5 GeV. The comparison with the 
^ recent experimental data is presented. 
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. I. INTRODUCTION 

i-C ' In the local NJL model [H, 0, [H meson interactions are described in terms of quark loops. If one neglects the 
dependence of external momenta in the corresponding integrals, the result will not violate chiral symmetry. In this way 
Q_ | i we can reproduce effective chiral Lagrangian corresponding to U(3) x U(3) symmetry Q, Q. In this lagrangian 
D | strong interaction vertices are expressed in terms of logarithmically divergent parts of corresponding loop amplitudes. 
Radiative interactions of mesons are described in terms of quark Feynman diagrams of anomalous type which do not 
contain ultraviolet divergencies. We use these anomalous quark loops for description of radiative mesons decays. 
t^J- . However for production of pseudoscalar and vector mesons in electron-positron collisions we should keep the external 
£>■ momentum dependence of the amplitudes. In this case, we encounter a serious problem of providing quark confinement 
condition. To solve this problem, one usually uses a nonlocal version of the NJL model, which involves the relevant 
formfactors for description of the interaction between mesons and quarks. The choice of formfactors leads to the 
functional unambiguity. Among the different ways to introduce these formfactors we should mention the QCD 
approach Q and the formfactors which arises from the instanton model Q. It is necessary to mention the models 
suggested by G. Efimov Q, Yu. Simonov @, Roberts and so on. 

It should be noted that a different form of formfactors leads to a different behavior of amplitudes in the physical 
region. This becomes essential at large values of external momenta. It is the reason why we apply rather a simple 
and rough method which consists in exact calculation of amplitudes and neglection of their imaginary parts to avoid 
the production of free quarks. As a result, we obtain a rather satisfactory description of radiative decay widths 
of vector and pseudoscalar mesons. This fact allows us to hope that this approach can be applied to describe the 
processes ee — » rf p, ee — > rfrc + Tr~ which can be measured in a series of existing and planned experiments with colliding 
electron-positron beams [l(|, [ll[, [T^, [T3 |. 

For the problems of the last type we hope to obtain only qualitative results in the center of mass energy range 
1—3 GeV. In our approach we do not introduce any quark-meson formfactors. 

Describing the decays of light mesons we ignore the dependence of corresponding loop integrals of external momenta 
- Approximation I. The problem of confinement in that approximation is automatically solved. In the case of heavy 
mesons we keep exact external momentum dependence of relevant loop amplitudes and neglect a possible imaginary 
part - Approximation II. Approximation II is used further for description of processes at electron-positron colliders. 
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II. RADIATIVE DECAYS OF VECTOR AND PSEUDOSCALAR MESONS 

For the description of interaction of mesons with quarks we use the NJL model lagrangian 0, H[ : 



eQA+ (ij 5 ) (g u \ u r] u + g s \ s rj s ) + ^ (^\ 3 p Q + \ u u + \ s £j q, (1) 



where q = (u,d,s) where u, d, s are the quark fields, Q — diag (2/3, — 1/3, — 1/3) is the quark charge matrix, 
\ u = (\/2Ao + As) /\/3, A s = (— A + v^Ag) /\/3 where Aj are the Gell-Mann matrices and A = a/2/3 diag(l, 1, 1). 
g u = m u / fir, g s = m s / ' f s arc the meson-quark coupling constants which are evaluated by Goldberger-Treiman relation 
(m u — 263 MeV, m s = 407 MeV are quark masses [141 ] ■ and /„• = 92.4 MeV is the pion decay constant and f s = 1.3/ ff ). 
g p = 5.94 is the p — > 2n coupling constant. 

Physical states of n and 7/ mesons are obtained after taking into account of singlet-octet mixing of rj u and rj s with 
the angle 6» = 51.3° 

■q = sin 6> + cos6>, (2) 
v' = Vu cos 9 + rj s sin 9. 

(3) 

We will consider the following processes below: p{uj) — > 777, 77' — > p(w)7, — ► viv')"/- 
The vector meson decay 

f (Pi) - 7(P2) + ^(Ps), 

is described by the amplitude of one loop with quark (see Fig. [TJ: 

M v ^p 1 = —^-2 e g P g v C PV M q J{p 2 1 ,0 7 pl)(e 1 e 2 piP2), (4) 
(2tt) 

where (abed) = e a /3 1 sa a b l3 c y d 5 , <jv = 5p/2, gp = g u if light quarks go through the loop and gp — g s if strange 
quarks are involved; Cpy is the flavour-color multiplier corresponding to quark- meson interaction, C nu — 2sin#, 
C np = 6 sin 9, C n ' U = 2cos0, C V > P = 6 cos 9, — 4 cos 9, C^^ = 4sin#; M q is the loop quark mass and 



J(pIpIpD = Rc (/ 



dk 1 



iir 2 {Ml -k 2 - i0) (Af g 2 - (fc + p 2 ) 2 - »0) (A/ 2 — (k — p 3 ) 2 - tO) 



= Re ( / da; / e?y — - 5 5 5 J , (5) 

\./o JO M g - ~ V Z P2 ~ ^3 ~ z0 / 

where z = 1 — x — y. In the heavy quark approximation (Approximation I) we obtain 

•HpIpIpD = (6) 
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Decay 


Experiment 


Approximation I 


Approximation II ((5]l 


p -*• VI 


39.47 


65 


33.72 


lu — > 777 


4.07 


7.83 


4.16 


7?' p7 


59.68 


76.18 


41.09 


77' — > LJ7 


6.15 


7.59 


4.04 


(/> -> 777 


55.59 


71.01 


117.9 


<^> -> ?7'7 


0.265 


0.497 


0.294 



TABLE I: The table of radiative decays. The values of the widths are in KeV. Approximation I - neglect of the external 
momentum dependence. Approximation II - the real part of exact loop integrals is taken into account. 

In this approximation a wide set of decays of light mesons was described in Q and the results were found to be in 
good agreement with the experimental data. The matrix element square can be written in the form: 



\M- 



V^P-y 



2 _ e 2 g 2 P g 2 v C PV _ ^ j^z^^y 



The phase volume of the final state is: 

And then the decay width reads as: 

1 a 



d$p 7 = 



v - M 2 P ) 



fp2fp3 _ \_M\ 
IE2 2i?3 



Ml 



8tt M 2 



V^P-f 



3 2 7 tt 4 V My 

The relevant expression for radiative pseudoscalar meson decays P — > V7 has the form 



a (M 2 ~M 2 ^' 



2 7 tt 4 V M P 



[gp 9v C PV M q J(M 2 ,0,M 2 )Y 



(7) 



(8) 



(9) 



(10) 



In Table Q] we present the theoretical results for both the methods - Approximation I ^ and Approximation II (J5]) 
- and compare them with the relevant experimental data. 

In particular, we would like to note that the ratio Rth. — ^{<t> — > v'l)/F(4> - * Wl) — 2-49 x 10^ 3 is in qualitative 
agreement with the result of the KLOE collaboration R exp . = (4.70 ± 0.47(stat.) ± 0.31 (syst.)) x 10~ 3 [13]. 

III. ASSOCIATIVE PRODUCTION OF PSEUDOSCALAR AND VECTOR MESONS IN 

ELECTRON-POSITRON ANNIHILATION 



The matrix elements of the processes of associative production of pseudoscalar and vector mesons 

e+(p + )+e-(p_)^V( Pl )+P(p 3 ), (11) 

where s = (p + +p_) 2 , p\ — m 2 , p\ = My, p 2 — Alp, in the lowest order of the QED coupling constant a have the 
form (see Fig. [5]): 



M PV = M*jmjA» 
S M 



where the QED lepton current is J^ m = v(p+)"/^u(p~) and the anomalous current has the form 

jA , = 9p gv Cpv (wip2) Mq J(p?)S)p 2 }) 
(2tt) 



(12) 



(13) 



where p+ + p~ = P2 = Pi + P3 and ei is the polarization vector of the vector meson (i.e. (eipi) = 0). 
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FIG. 2: The processes at electron-positron colliders. 



The cross section built by general rules is 



da 



— Vim 

8s 



PV, 



where the phase volume of the final state has the form: 



PV 



d 3 pi d 3 p 3 1 
2#7 2E 3 (2tt) ; 



-<5 4 (p+ +p- -p! -p 3 ). 



(14) 



(15) 



As we are concerned with the total cross section only we can use the property of anomalous current gauge invariance 
and thus rewrite the final state phase integral as 



(16) 



The second term in the braces does not give a contribution due to gauge invariance of the lepton current J^ m . The 
first term contribution is proportional to 



^2 (einpip 2 ) (e*^pip 2 ) = ~X (s, Mp, M, 



(17) 



where A(et, b, c) 



2 . 



b +cr — lab — 2ac~ 26c is the well-known triangle function. Thus, the quantity | Jrf | in |T6 



does not depend on the vectors p\ and P2 themselves but only of their squares p\ = My, p\ = Mp. This allows us to 
calculate the integral over the final state phase volume which, neglecting the masses of leptons, can be written as 



d 3 pi d 3 p 3 



2Ei 2£^3 

Then the total cross section obtains the form 

.,2 



$ (<1 ~ Pi ~ Pa) = 7; 



7t (s,Mp,M v ) 



<r(s) = ^^aI (a, Mp, My) \g v g P C PV M q J (M P , My, s)\ 2 . 



The differential cross section can be written as: 

da 

lm 2 



a(s) 



3(1 + cos 2 0) 



16tt 



(18) 



(19) 



(20) 



where 9 is the center of mass angle between the direction of 3-momenta of the initial electron p- and the final vector 
particle momentum direction p\ . 

Let us consider the concrete process ee — > r/p. The expression for the total cross section (fT9|) implied only contact 
Feynman diagram, i.e., Fig [5] a. Recalling the possible conversion of virtual photon into vector mesons beyond 
resonances one must take into account the diagrams presented on Fig[^l b. This leads to the replacement of the factor 
s~ 3 in p^|) by the following one: 



1 



1 

2(1- 



(21) 
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FIG. 3: The comparison of our result for the e + e — > r]'ir + -K with the BABAR-collaboration results for the e + e — > r/7r + 7r 
channel Hill. 



The cross section of the process ee — > r/(950)p is drawn on Fig. [31 where the relevant experimental data are also 
shown. One can conclude that satisfactory agreement within the experimental errors is observed. 
Let us now make a prediction for the process ee — > rf<p. The relevant correction factor is 



1 



/ 



V 



1 



(22) 



Besides only s-quark loop works. The result is given in Fig. fH 



IV. CONCLUSION 



In this paper, we investigated the radiative decays of vector and pseudoscalar mesons described by the quark loops 
of anomalous type. 

Let us note that in [l6| we considered the process <j) — > /q (980)7 within the same framework of the NJL model. 
However, there the quark loop contribution was small enough and the main contribution arose from terms of next 
order of l/N c expansion (where N c is the number of colors) - meson loops. In this paper, we have another situation: 
meson loops are absent totaly and only quark loops of anomalous type give a contribution to the amplitude of the 
process. 

Both the approaches (Approximation I and Approximation II) were considered. We show that the application of 
the NJL model leads to rather satisfactory agreement with the modern experimental data for the radiative decays. 
That allows us to use the Approximation II to calculate the cross sections of associative vector and pseudoscalar 
mesons production in the electron-positron annihilation channel in the lowest order of electromagnetic constant. In 
the last case, the heavy virtual photon converts to the set of pseudoscalar and vector mesons. Two mechanisms 
must be taken into account: first one with the intermediate virtual photon and the second one which contains the 
conversion of intermediate photon into vector meson. We give a comparison of our prediction for the process ee — > r/' p 
with the experimental data of the BABAR collaboration. For the process ee — > rf<p the prediction is given for future 
experimental data. 
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FIG. 4: The prediction for cross section e + e 



For comparison with the experiment the precision of our results is worth mentioning. We should like to notice that 
for the NJL model results precision is of an order of 20-30 %. 

One of our important theoretical assumptions is the absence of the imaginary part of relevant amplitudes. The 
mechanism of elimination of the imaginary part is tightly connected with the confinement nature and is not considered 
here. We carry out the elimination "by hand" ("naive confinement"). It is to be noted, however, that if the imaginary 
part is taken into account, the considerable disagreement with the experimental data will occur in decay case. For 
instance, if the imaginary part of the amplitude is taken into account, the decay width of </> 



7/7 is T th 



0.824 KeV while the experiment gives T 



exp. 



*V 7 



0.265 KeV (see Table QJ. 



Concerning the singlet-octet angle mixing we use the additional interaction of the t'Hooft type in lagrangian in the 
NJL model [l5j] , [22j , [HI . This approach was widely used in literature (24| , (25| . 

Let us note however that the alternative solution of mixing angle problem was developed in [lj| , [H| , [13] , [HI • 
In this approach two mixing angles appear. The application of the both approach to describe the decays of the 
pseudoscalar and vector mesons leads to the similar results. 

Our results for decays in Approximation I are in agreement with the ones obtained in [l7| (compare Table [J and 
Table 2 in [13). 
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